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Abstract 

I give an efficient algorithm for the reduction of multi-leg one- loop integrals of rank one. 
The method combines the basic ideas of the spinor algebra approach with the dual vector 
approach and is applicable to box integrals or higher point integrals, if at least one external 
leg is massless. This method does not introduce Gram determinants in the denominator. 
It completes an algorithm recently given by R. Pittau. 
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1 Introduction 



Next-to- leading order calculations have become nowadays the standard for perturbative calcu- 
lations in high energy physics. One of the major obstacles is the tedious reduction of tensor 
loop integrals (e.g. integrals where the loop momentum appears in the numerator) to standard 
scalar integrals. General algorithm have been known for a long time (like the Passarino-Veltman 
algorithm or the Feynman parameter space technique ||). Although they solve the prob- 
lem in principle, the resulting expressions are usually quite long and involve the appearence of 
unphysical singularities in the form of Gram determinants in the denominator. These singular- 
ities usually drop out only after combination of many different structured terms. 
A variety of alternative approaches exists: J.M. Campbell, E.W.N. Glover and D.J. Miller || 
introduced scalar integral functions in higher dimensions and absorbed the Gram determinants 
in a linear combination of these. D.B. Melrose [|J and independently W.L. van Neerven and 
J. A.M. Vermaseren || have used the fact, that in four dimensions there are maximally four 
linearly independent vectors, to derive a relation between a scalar pentagon integral and scalar 
box integrals. The derivation of van Neerven and Vermaseren was based on dual vectors. The 
dual vector approach was further used by G.J. van Oldenborgh and J. A.M. Vermaseren || and 
A. Signer ||. Recently R. Pittau |7j has given an algorithm, which makes use of the spinor 
algebra, to reduce tensor integrals to scalar and rank one integrals. In this paper I give a 
formula which allows the reduction of the remaining rank one integrals to scalar integrals. It 
combines the basic ideas of the spinor algebra approach and the dual vector approach. This 
formula completes the algorithm given by R. Pittau. 

This paper is organized as follows: The following section introduces the necesarry notation. In 
section 3 I derive the basic formula. The complete algorithm is given in section 4. An example 
is worked out in section 5. Section 6 contains the conclusions. 



2 Notation 



For an introduction to spinor products I refer to [11]. A typical one- loop integral has the form 

d D k ^ ± \..4...\ qi ±)...(p m ± \..4...\q m ±) 



(1) 



where the Pi and the are external massless momenta and the iV, denote the internal prop- 
agators. The loop momentum k is allowed to appear more than once in each sandwich 
(Pi ± |...$...|(&±). This form is general: After the decay of (possible) intermediate heavy parti- 
cles all external particles can be taken to be gluons, photons or light fermions. The polarization 
vectors of gluons and photons are expressed in the above form using the spinor helicity method 
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where q is an arbitrary null reference momentum. It will be usefull to introduce a slight 
generalization of spinor products: Consider first elements of the form ^ = a/^. They form an 
algebra, which I will denote by A. It is convenient to think of the space of ket-spinors generated 
by the spinors obtained from all external momenta as an left ^-module V. The dual space of 
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bra-spinors is denoted by V* and is a right A-module. The spaces V and V* decompose into a 
direct sum of subspaces of positive and negative helicity: 

V = v+ © y_ V* = V* © K! (3) 

Non-vanishing spinor products are only obtained between elements of V? and or between 
elements of V* and V+. They are denoted by 

[06] = (a + |6-) (06) = (a - |6+) (4) 

where (a ± | G V± and |6±) G V±. Multiplication of a spinor with an element of A changes the 
helicity, e.g. if \p+) G V + and </l G ^4 then G It is convenient to define a "degree" of 

these objects. Assume that the spinor \q) corresponds to the external momentum q and that 

\p) = (/i n ..4i\q). (5) 

Then we define deg(|g)) = and 

degflp)) = n. (6) 

A similar definition applies to bra-spinors. Changing the helicity or flipping a spinor from bra 
to ket (or vice versa) does not change the degree. It follows that for arbitrary spinors \p), \q), 
\k) and \j) 

(qp) = 

[qp] = -(-l)^P(-l)^^\pq], 
(pq)(kj) = (pj)(kq) + (-l)^%-l) d ^ k {pk)(qj), 
[pq}[kj] = \pj}[kq} + (-^ de9q (-^ 9k \pk}[Q3}- (7) 
For deg p = deg q = deg k = these formulae coincide with the usual ones for spinor products. 

I take the sign of the antisymmetric tensor as £0123 — +1- It is useful to introduce the fol- 
lowing short-hand notation: s(q a , q b , q c , q d ) = ^ie^ pa q^qlq p c q^ 



3 Derivation of the basic formula 

Combining the basic ideas of reduction methods based on dual vectors and on the spinor algebra, 
I derive a formula, which allows an efficient reduction of higher point integrals. 
The kinematics are shown in fig.0 The pi are supposed to be external momenta whereas k is 
taken as a loop momentum. I assume that the pi are linearly independent. This means that 
the method applies to pentagon or higher point integrals. A similar formula for box integrals 
is given later on. Let us introduce the notation 

i 

J'=l 

ki = k- q { . (8) 

It is convenient to set k = k and ^5 = — q^, so that one basically always deals with five external 
momenta pi, ...,p 5 , even for hexagon or higher point integrals. We have 

2k ■ ft = k 2 - kf + qf, 

2k - Pi = kU-kt + <fc 2 -<&.!■ (9) 
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Figure 1: A subgraph of a one- loop diagram. The loop momentum is flowing through k and 
fc 4 , the p^s are assumed to be the external momenta. 



It is convenient to work in a scheme where the Dirac matrices are in four dimensions (like the 
four-dimensional helicity scheme [Q). If the integral is finite we can take all vectors in four 
dimensions. The case of a divergent integral needs some special care: I follow the approach 
by G. Mahlon Jl3| and split the loop momentum into a four-dimensional part and a (— 2e)- 
dimensional part 

fcP) = fc(4)+jfe(-*). (10) 

It is further assumed that a four-dimensional vector is orthogonal to a (— 2e)-dimensional one. 
I start from the dual vectors written as 

< = \e ahcd ^e^q v M, (11) 

and the Schouten identity 

e(qi,q2,q3,q4)(k^f = £ (k ■ v a ) (k ■ q a ) . (12) 

a=l 

This equation is related to four dimensions, therefore k^ appears on the left hand side. Since 
q a and v a are four-dimensional objects, the scalar products of k^' 2 ^ with q a or v a vanish and I 
dropped therefore the superscripts. Now 

k-v a = ±\e abcd {2Ti ± U h U d -T4i4 c i d ) 
o 

= ± l -e abcd 2Tr ± MMd (13) 
o 

where Tr± denotes a trace of Dirac matrices with the insertion of a projection operator 1 (1 ± 75). 
Since 

Tr{Ubicid) = {2k ■ q b ) {2q c ■ q d ) - (2k ■ q c ) (2q b ■ q d ) + (2k ■ q d ) (2q b ■ q c ) (14) 

and each terms has one scalar product symmetric in two indices, these terms drop out when 
contracted with the antisymmetric tensor. Working out the r.h.s of eq. ([12]) yields 



(k ■ v a ) (k -q a ) = ± (2k ■ q{Ti ± M2fah ~ ^k ■ q{^MM± 

+2k ■ g 3 Tr±^i^4 - 2k ■ q^±Uihh) ( 15 ) 
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Writing 2k ■ q a = k 2 — k 2 + q 2 , one sees that the first two terms cancel propagators, whereas 
the last term yields a contribution of 

= ± {Ti ± Uihhhh + qlqlikl - kl) + q \q\(kl - k\) + q\q\{k\ - k 2 )) (16) 

The r.h.s. of eq. ( |i6|) has been arranged such that no constant term survives. One finally 
obtains 

-h?Tr ± 1kofafafa - {k^f (Tr ± A^4 - Tr T ^^ 4 ) • (17) 

This equation holds if the Dirac matrices and the external momenta are in four dimensions. The 
loop momentum k may be in D dimensions, allowing the application also to divergent integrals. 

A careful inspection of the basic equation ( |Ti2| ) shows that the term e(qi, q 2 , (fe, q±) on the 

left-hand side has to be in four dimensions, whereas the manipulations on the right-hand side 
involved only scalarproducts 2k ■ q a = k 2 — k 2 + q 2 which may be continued into D dimensions 
as long as the external momenta q a stay in four dimensions. The last term of eq. fllTD is a 
correction term, which takes care of divergent integrals. If the integral is finite, k may be taken 
to be four dimensional and the last term vanishes. It will turn out that the correction term is 
irrelevant in the application of eq. fll7D to pentagon or higher point integrals. The equation 
above may be rewritten as : 

IM#M&M) = "Tr T ijkfcfafa) (No + M 2 ) 

--(B ±e(k- Pl ,p 2 , P3 , Pi )) (N Q + M 2 ) 

+^ (B 1 ±e(k, Pl + P2 ,p 3 ,p 4 )) (iVi + M 2 ) 

-^(B 2 ±e{k, Pl ,p 2 +p 3 ,p 4 )) (N 2 + M 2 2 ) 

+- (B 3 ± e(k, Pl ,p 2 ,p 3 +pt)) (N 3 + Ml) 

-i (B 4 ± e(k, Pl ,p 2 ,p 3 )) (iV 4 + M 4 2 ) 

(18) 

with Nt = kf — M 2 , where the internal propagators and masses are denoted by AT, and Mj, 
respectively. Bi depends only on the pinched integral under consideration and is given by 

Bi = st — m 1 m 3 — m 2 m 4 . (19) 

Starting from the pentagon integral the internal propagator % is first pinched. This gives a box 
integral and the value of Bi is calculated from that box integral according to eq. ( JT9| ) and 
according to the kinematics shown in fig. |2|. Clearly the value of B{ does not depend on the 
starting point for the labelling of the external legs lj of the box integral. Formula ( |I9"D also 
applies to hexagon or higher point integrals, where the sum of the external momenta from leg 
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/ \ s = (h + l 2 f t=(l 2 + h) 2 



Figure 2: The labelling of a box integral with external momenta Zj. 




Figure 3: The labelling of a triangle integral. 



5 to leg n is identified with p 5 . 

For a tensor box-integral one can use 

-\c G (No + Ml) + ±a (N, + Ml) 

-\c 2 (N 2 + Ml) + itf 3 (N 3 + Mi) (20) 

with 

Q = m\ + ml- m\. (21) 

This equation is just obtained by expanding the trace. The orientation of the triangles de- 
pends on how the triangles are obtained from the box integral : If the loop momentum k Q 
is not pinched, j = ko and ml is the first external mass in the direction of the flow of j. If 
fco is the pinched momentum, then j = k\. The kinematics are also shown in fig|| Equa- 
tion (|IBD and (|20| ) are the basic equations for the reduction of tensor loop integrals. For box 
integrals the terms with the antisymmetric tensor in the numerator will vanish after integration. 

It remains to investigate under which condition the traces may be formed. It turns out that it 
is always possible to form a trace Tr ± ^!^2^3^4^5 with p 5 = —pi — p 2 — p 3 — p 4 , if at least one 
external leg is massless. Terms like (q — ), where k is the loop momentum, appear in the 
numerator of tensor integrals . Using the Schouten identity twice, one can write: 

[ij\(ij)(q-mi-) = -(-i) de9i (-i) de9j (d)\jw-mi-) 
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-{-I) de9i {-l) de9j {qi)\il]{3-m3-) 
+(qj)[il)(i-mj-) 

+(-l) de9i (-l) de9j (qi)[jl](i + mj+). (22) 

This will be the basic equation (together with suitable choices of % and j) for forming the traces. 
If a n-point integral (n > 4) has at least one massless external leg, say p%, one can choose 

= \Pi~), 

= hhhh\pi-) (23) 

where p$ = —p^ — P3 — p 2 — Pi- Of course, if there are more massless legs, simpler choices are 
possible. It is easily checked by induction that the terms (i — and (j — always 

cancel propagators. 



4 The algorithm 



This gives us the following algorithm for the reduction of multi-leg loop integrals (e.g. at least 
four external legs), which is applicable if at least one external leg is massless. Three-point and 
two-point integrals are usually not so complicated and can be reduced using standard tech- 
niques. 

Step 1 consists in reducing tensor integrals to rank one integrals and was given by R. Pit- 
tau 0: If there at least two massless legs, say p% and p 2 , one writes 

$ = 7 ^—[(2k-p 2 )p' 1 + (2k-p 1 )^-^ 2 -^ 1 }. (24) 

*Pl " P2 

The first two terms already cancel propagators, whereas a product of the last two terms can 
be rewritten as 

(l-|#|2-)(2-|fll-) = (2k-p 1 )(2k-p 2 )-(2p 1 -p 2 )k 2 , 
(l-|fl2-><2 + |Jf|l+) = 1 [(2 Pl -p 3 ) (2 ft -A;) (l-\m-) 

- (2 Pl ■ k) (i - 

+ (2p 3 -A;)(l-|^#i|2-) 

+^(1-1^3^12-)], (25) 

where p% is another (not necessarily massless) external leg. If there are no massless legs, or only 
one massless leg, one can take linear combinations of the external momenta in order to form 
two massless momenta [|]]. If all external momenta are massive this will re-introduce Gram 
determinants. 



In step 2 we reduce the remaining rank one integrals to scalar integrals by first forming the 
trace with the help of equation (|22). This can always be done if at least one external leg is 
massless. Using equation (]18D and equation ( f20|) we then obtain scalar integrals, integrals with 
an antisymmetric tensor in the numerator and a correction term. For box integrals the terms 
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Figure 4: The Feynman diagram with the helicity assignment. All external momenta are taken 
to be outgoing. 



with the antisymmetric tensor in the numerator will vanish after integration. For higher point 
integrals one may rewrite 

e(k,p h p2,p 3 ) = Tr + (^i^ 3 ) - Tr_ (#1^3) (26) 
and iterate the procedure. 

Step 3 : The correction term is evaluated as follows [ 13 ]: 

T^Fv xT . (27) 



(27r) 4 -^iV 1 Ar 2 ...Ar n v 7 7 (27r) 6 - 2£ ^ 2 ...A^ n ' 

The correction term is proportional to the scalar integral in D = 6 — 2e dimensions. These 
integrals are finite in the limit e — > for n > 4 || and the correction term does therefore not 
contribute to order 0(e°). 



5 An example 

As an example I consider here a pentagon integral which was recently encountered in the one- 



loop calculation for e + e — > qqQQ fl5[. We want to evaluate the helicity diagram shown in fig. 



g Thedia gram corresponds to the following integral: 



[12] (34) r d D k (5-IW0I6-) , . 

/ OttW k 2 k 2 k 2 k 2 k 2 ^ ' 



s 56 J (2tt)^ fc 2 ^!^ 
The numerator is rewritten as 

(5-]JMMM6-> = -(5-1(^3+^)^(^+^)16-) 

+k 2 2 (5 - \to\6-) - k 2 2 (5 - (29) 

The first term on the r.h.s of eq. fl29|) corresponds to a rank one pentagon integral, the second 
term to a rank one box integral, whereas the last term already corresponds to a scalar box 
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integral. The rank one pentagon integral is reduced as follows: Multiplying with [23] (23) and 
using eq. ( p2]) we obtain: 



[23] (23) (5 - \(fa + fa)h(fa + fa)\Q-) = 

= -(6 + |1 + 2|3+) [34] (45) (k\ - kfj - [61] (12)(2 + |3 + 4|5+) (k 2 2 - k 
(6+ |l + 2|3+)(2+ |3 + 4|5+) 



(34) (4+ |5 + 6|1+) [12] 
[61] (45) 



-Tr + (f 2 fafa (fa + fa) fa fa 



(4- |5 + 6|1-) 
According to eq. (|i~8D the traces give after integration 

d D k Tr±(hfafa(fa+fa)fafa) f d D k 



^-(hfafa(fa+fa)fafa)- (30) 



(9ir} D h 2 h 2 h 2 h 2 h 2 / Ott) d h 2 h 2 h 2 h 2 h 2 

y£t l\ J A, q i\ A, 2 A, ^ J \ ) 12 3 4 



-k 2 2 i:i T (fa fa (^5 +fa)fa 



— 2 ( s 123 s 234 — •523 s 56) ^2 + T^ 12 " 5234 ^ 3 ~~ T^ 12 " 523 ^ 4 + ^"^S^ 34 ^) — T^ 123 " 534 ^! 

+0(5) (31) 
The rank one box integral is evaluated as follow: After multiplication with [14] (14) we obtain 

[14]<14><5- |JS |6-> = 
= -(54) [46] (k 2 Q - k 2 ) - (51) [16] (kl - k\ + s 56 - si23 

Tr+ (Ui(J>2+fa)fa) 



(4- |2 + 3|1-) 
(51) [46] 



Ti-(hfa(fa+fa)fa) (32) 



(4+ |2 + 3|1+) 
The traces give after integration (using eq. Q2l]|)): 

I d D k Ti±(hfa(fa + fa)fa) = i r d D k i (si23S234 

- («23 - S234) ^0 + ( S 56 - S123) ^1 

- (^56 - S234) ^3 + ( S 23 - S123) ^4) (33) 



This completes the reduction to scalar integrals. In this example all propagators were massless. 
The extension to massive lines is straightforward. If for example &2 is replaced by a massive 
line, the algorithm yields a combination of scalar boxes and a scalar pentagon integral. The 
scalar pentagon integral can be expressed in terms of scalar boxes according to well-known 
formulae P, HI. 



6 Conclusions 

In this paper, I have given a formula which allows the reduction of rank one pentagon or higher 
point integrals. If there is at least one massless external leg, this formula can always be applied. 
This method does not introduce a Gram determinant in the denominator. A corresponding for- 
mula for rank one box-integrals was also given. For rank one triangles or bubbles conventional 
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techniques are usually efficient enough. The formulae in this paper complete an algorithm 
recently given by R. Pittau, and form together an efficient algorithm for reducing tensor loop 
integrals. 



I would like to thank D.A. Kosower and J. A.M. Vermaseren for useful discussions and for 
suggestions concerning the manuscript. 
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